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Abstract 

A cranked mean-field model with two-body T=l and T=0 pairing interac- 
tions is presented. Approximate projection onto good particle-number is en- 
forced via an extended Lipkin-Nogami scheme. Our calculations suggest the 
simultaneous presence of both T=0 and T=l pairing modes in N=Z nuclei. 
The transitions between different pairing phases are discussed as a function 
of neutron/proton excess, T z , and rotational frequency, Twj. The additional 
binding energy due to the T=0 np-pairing correlations, is suggested as a pos- 
sible microscopic explanation of the Wigner energy term in even-even nuclei. 
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The study of pairing correlations is one of a central theme in nuclear structure physics. 
Although the energy gain due to pairing correlations is rather modest, these correlations 
strongly influence many properties of the atomic nucleus. The large body of phenomena re- 
lated to pairing among like-particles can be well understood, at least qualitatively, in terms 
of the simple BCS model with seniority force. In light nuclei, especially those with N=Z, it is 
well established empirically that neutron-proton (rap) short range (pairing) correlations are 
of importance. The mean-field formalism like the Hartree-Fock-Bogolyubov (HFB) method 
is in principle capable to simultaneously treat both T=0 and T=l pairing modes. The nec- 
essary generalizations of the Bogolyubov transformation were worked out by Goswami and 
coworkers ]1|^3| and Goodman [|]]. These early calculations indicated that the importance 
of rap-pairing is restricted to the vicinity of the N=Z line, ||, and that the T=0 and T=l 
pairing phases are exclusive, see also |6|,[7j . The formalism was further extended to describe 
rotating nuclei || , suggesting a possible phase transition between the two pairing modes at 
high spin ||. 

The recent progress in nuclear spectroscopy, related to the event of highly efficient detec- 
tor arrays and development of radioactive ion beam facilities, is opening up new avenues to 
study the nature of nuclear interactions, in particular, rap-pairing correlations at the N=Z 
line. Phenomena like possible phase transition between different pairing modes in rapidly 
rotating nuclei, the influence of rap-pairing on the position of the proton drip line and the 
stability of drip line nuclei due to the additional binding energy emerging from rap-pairing 
are becoming important issues in nuclear structure. 

The aim of this paper is to investigate basic features of rap-pairing. We present a model 
applicable for pairing-and-deformation self-consistent cranking calculations and introduce a 
method to restore approximately the particle-number symmetry. This concept is an exten- 
sion of the so called Lipkin-Nogami method for the case of a non-separable proton-neutron 
system. The method is independent of the kind of two-body interaction used in the calcu- 
lations. Applying approximate number-projection, results in the simultaneous presence of 
both T=0 and T=l pairing modes. A detailed discussion of different aspects of our model 
will be given in a subsequent publication. 

The starting point of our calculations are the eigenstates of a deformed phenomenolog- 
ical single-particle potential and therefore, spherical symmetry is broken already from the 
beginning. The main goal is to construct a formalism which is flexible enough to account 
for simultaneous scattering of (?) the nucleonic pairs where both particles occupy states of 
different signatures and (ii) pairs where both particles occupy states of the same signature. 
According to Ref. ||, these two pairing modes will be further denoted as a/3 and ct(3, 
respectively. 

The most general Bogolyubov transformation can be written as: 
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where a (a) denote single particle states (including isospin indices) of signature r = 
respectively, while j labels quasiparticles. Following the calculations of Ref. 0, we will 
further impose the so called antilinear simplex symmetry, = PTR Z as a self-consistent 
symmetry (SCS), see also One should bear in mind that, due to the antilinearity of S A 



the transformation properties of creation and destruction operators with respect to Sf will 
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depend on the phases of the single-particle states. In other words one cannot introduce any 
new quantum number associated with that symmetry. After applying the S£ symmetry the 
Bogolyubov transformation still remains complex but the imaginary part decouples from the 
real part in the sense that the different signature blocks of the density matrix, p = V*V T , 
and pairing tensor n = V*U T are either real or imaginary: 
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Furthermore, the complex structure of the single particle potential, h, and the pairing 
potential, A, and consequently the HFB equations are fully determined by the complex 
structure of the p and k matrices, respectively. 

In this work we restrict the two-body np-pairing interaction to a simple extension of the 
standard seniority pairing interaction. It is separable in the particle-particle channel, v a/ 3 y s oc 
g a pg*$, with g a p proportional (up to a phase factor) to the overlap (a T \/3 T r) between single- 
particle wave functions. Apart from weak modifications due to the isovector components of 
the nuclear-mean field, like the static Coulomb potential, the interaction is dominated by 
ad and aa types of pairing. 

In order to relate the structure of the act and aa pairing modes to the isospin quantum 
numbers, let us briefly consider a nucleus with isospin- and time-reversal symmetry. By 
decomposing the pairing potential into the different isospin components (T,T 2 ) one finds 
that the T=l and T=0 components of the np-pairing depend on the combinations of the 
same elements of the pairing tensor but with opposite sign [ffl. Consequently, with the 
pairing tensor of the form of (0), the T=0 component of aa pairing is ruled out due to the 
Sf symmetry. Similar analysis shows that T=l component of the aa pairing also vanishes. 
Therefore, in our model the act pairing is equivalent to T=l and aa to T=0. This simple 
analysis also reveals the important role played by symmetries in the theoretical description 
of np-pairing, see also ||. 

The average gap parameters are equal to (we adopt the convention where r = 1 for 
neutrons and r = — 1 for protons) 



A:;i = -^A^ where = G'^ £ K(Xr ~, (4) 

a T >0 



for T=l pp- (nn-) pairing, 
A I:fc = -<«^-> A np =1 where A^ =1 = V =1 £ (a n \(3 p ) [k^ + Kpp ~, } (5) 



a n ,/3 P >0 



for T=l (aa) np-pairing and 

A£3_ T = tr(a T \(3^(A T n ^) and Ag*L- = -.r(a T |/3_ T )3(A^= ) 
where A£ p =° = |G£=° £ K|/3 p ) {^(k-^ - S^j} 

a n ,(3 p >0 
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for T=0 (aa) rap-pairing. The strengths of the interaction are denoted by G> T /. 

To prevent a sudden collapse of the static pairing correlations, e.g. induced by fast 
nuclear rotation we introduce an approximate particle-number projection using the Lipkin- 



Nogami (LN) method |IT |. This method is equivalent to a restricted HFB-type variation, 



5{HFB\H"\HFB) = 0, for the Routhian: 

H" = H U} -J2 \^AN T - ]T \®AN T AN T/ . (7) 

In the LN method the parameters A^ 1 -* are standard Lagrange-type multipliers whereas the 

(2) 

parameters A^ T / are kept constant during the variational procedure and eventually adjusted 
self-consistently using three additional subsidiary conditions. 

(H"(AN T AN T , - (AN T AN T f))) = 0. (8) 

where AN T = N T —N T and the symbol ( ) stands for the average over the \HFB) state. The 
LN theory is technically similar to the HFB theory but for the Routhian (0). The resulting 
LN equations take the form of HFB equations with a single-particle field and pairing field 
renormalized as follows: 

h% -> h TT , + 2ASW and -> A TT , - 2\®k tt ,. (9) 

An open question in mean-field calculations with rap-pairing is related to the strength 
of the interaction. Whereas the strength of the pp- and nn- seniority pairing force is well 
established by a fit to the odd-even mass differences, very little is known about the strength 
of the rap-pairing force. Based on isospin-symmetry arguments, it seems well justified to 
assume that at the N~Z line Gp^n) ~ G^p 1 - Therefore, our results will be presented either 
as a function of or at a given value of the parameter x T=0 that scales the strength of T=0 re- 
pairing with respect to the average strength calculated for nn- and pp- pairing correlations 
i.e. x T =° = Gl;VGl;' while G^ 1 = (G^ 1 + G^)/2. 

Fig. 1 shows the pairing gaps at zero frequency for a self-conjugate, N=Z, nucleus 
calculated with (BCSLN) and without (BCS) approximate number projection as a func- 
tion of x T=0 . In this case we also disregard the Coulomb interaction and, therefore, 
^pp =1 = ^un 1 = ^0 and Gpp 1 = G^ 1 = G^ 1 . The BCS version of our model has 
been discussed in the literature [|J and the solutions can be characterized as follows: (i) 
For x T=0 < 1 (G^p < G^p 1 ) the T=l pairing is energetically favoured over the T=0 pair- 
ing. The pairing energy depends only on A 2 = 2Aq + (A^ 1 ) 2 and no energy is gained 
by activating the T=l rap-pairing, (it) The solution at x T=0 = 1 (G^ = G^" 1 ) is highly 
degenerate. The HFB energy depends only on A 2 = 2Aq + (A^ 1 ) 2 + |A^=°| 2 . Also in this 
limit, no energy is gained due to rap-pairing, (in) The solution at x T=0 > 1 (G^ > G^" 1 ) 
corresponds to a pure T=0 rap-pairing phase. 

As shown in Fig. 1, the results from the number-projected calculations are quite different. 
The critical value of the strength necessary to activate T=0 rap-pairing is larger, x T c ~^ ~ 1.1. 
The LN- method introduces different modifications of the pairing potential for the T=l pp- 
(nn-) and T=l rip-pairing field. It causes that at x T=0 < x T c ~^ and under the assumption of 
Gnp 1 = G^r 1 the pairing gap for the T=l rrp-pairing, A^J 1 , becomes zero. It implies that 
the LN method has an isovector component, that requires further investigation. However, 
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at x T=0 > x^~^ the T=0 np-pairing correlations coexists with the T=l, |T Z | = 1 pairing. The 
exclusivness of T=0 and T=l pairing phases in N=Z nuclei is a generic feature of the BCS 
method and is smeared out in the number-projected calculations. 

For N^Z both BCS and BCSLN models provide solutions which are qualitatively similar 
to the T z =0, BCSLN case i.e. the T=0 np-pairing correlations coexist with the T=l pp- and 
nn-pairing correlations. It is worth stressing that the value of x^7i° ^ s st ron gly T 2 dependent 
i.e. increases quite rapidly with neutron/proton excess or, alternatively, np-correlations are 
restricted to small |T 2 |. In consequence, our calculations suggest a critical value of |T^"*| 
beyond which, i.e., for \T Z \ > |T° n *|, there is no collective solution to np-pairing, see Fig. 2. 

The additional binding arising from T=0, np-correlations as well as the narrow region 
where it is active, is shown in Fig. 2. Early calculations of nuclear masses based on a 
macroscopic-microscopic approach have shown particularly strong deviations from the ex- 
perimental data in the vicinity of the N~Z line [|l|]. Most probably, only part of these 
deviations can be attributed to np-pairing while part of it can be accounted for by the self- 
consistent mean-field. One can assume that the single-particle mean field is properly taken 
into account by the extended Thomas-Fermi model which is a semi-classical approximation 



to the Hartree-Fock method [13]. In this model there is a systematic binding energy offset 



by ~ 2 MeV at the N=Z line that one might associate to the lack of rap-pairing correlations. 
Based on this assumption and the results of [13] we can estimate the strength of the T=0 
pairing force G T=0 ~ 1.2G T=1 as shown in Fig. 2. The shell-model estimate of Ref. |14 



yields G T ~ 1.3G T 1 . Modern versions of macroscopic- microscopic mass calculations |T5[ 



cure the difficulties arising around the N~Z line by introducing an extra term to the liquid 
drop formula - the so called Wigner energy flEfl . Hence, a possible microscopic origin of the 
Wigner term in even-even nuclei are the T=0 np-pairing correlations. 

A third way to generate a phase transition from T=l to T=0 pairing is by rotation ||[|. 
The results of cranking calculation for 46 > 48 Cr are shown in Fig. 3. These are qualitative 
calculations at constant deformation and as such should not directly be compared to the 
experimental data. The calculations are performed with and without the T=0 force. At 
hu) = 0, the strength of the T=0 force is undercritical and T=0 pairing is not active. At a 
certain critical frequency, hu crit , there is a sudden onset of T=0 pairing, see also ||. The 
latter effect can be viewed as either a phase transition or band crossing. The coherent action 
of centrifugal and Coriolis forces tend to align the angular momentum of the quasi-particles 
along the rotational axis. It weakens the T=l pairing correlations and simultaneously, 
increases the number of pairs of nucleons with parallel coupled angular momenta, thus 
enforcing the T=0, ckck, np-pairing correlations. In the T=0 phase, angular momentum is 
built by smoothly aligning np-pairs along the rotational axis, without involving any pair 
breaking mechanism. This situation is totally different from the well known response of the 
T=l pairing to nuclear rotation, where pairs are broken to generate angular momentum. 
With increasing frequency, the T=0 pairing correlations tends to saturate. Note that the 
T=0 phase coexists with T=l pp- and nn-pairing phases although the onset of T=0 np- 
correlations suppresses the T=l phase. However, our calculations show that the T=0 np- 
pairing and T=l np-pairing phases are always exclusive at hu ^ 0, independently on their 
relative strengths. 

The influence of the T=0 np-pairing correlations on the dynamical moment of inertia is 
conspicuous. Nuclear rotation in the presence of T=0 pairing correlations resembles classical 
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rigid body like rotation even though T=l pp- and nn-pairing correlations are present. Note 
also that the moments of inertia at large frequency in the presence of T=0 pairing exceed by 
far the value obtained for a system without pairing, see Fig. 3. Even though the situations 
discussed above were visualized for 46 > 48 Cr only, these classes of solutions appear to be generic 
for all even-even nuclei with N~Z, i.e. do not depend qualitatively on A. 

Our results can be summarized as follow: The previously suggested exclusiveness of 
the T=0 and T=l pairing phases H0,fj does no ^ su PP°rt in our calculations. The 
sudden phase transition between the T=0 and T=l pairing modes is a generic feature of the 
BCS approximation for N=Z nuclei. This phase transition becomes smeared out in number- 
projected LN calculations. There, the T=0 np-pairing correlations coexist with T=l nn- and 
pp-pairing correlations over a broad range of the strength as well as rotational frequency 
hu, when x T=0 > x T c ~^ . However, pairing correlations of aa and act type counteract. For 
N^Z, the T=0 np-pairing correlations and T=l nn- and pp-pairing correlations do coexist 
in both BCS and number-projected LN calculations. The T=0 np-pairing correlations are 
confined to a narrow region along the N=Z line. The additional binding arising from these 
correlations may be viewed as a microscopic origin of the Wigner term in even-even nuclei. 
Even in the cases where np-pairing correlations are not present in the ground state one 
can generate T=0 np-pairing correlations at large rotational frequencies. An onset of aa 
np-pairing quite dramatically influences the mechanism of building angular momenta. The 
T=0 np-pairs can easily be decoupled from the deformed core by the Coriolis force and 
consequently nuclear rotation resembles rigid-body like rotation though pp- and nn-pairing 
correlations are present. In both BCS and number-projected calculations, the T=0 and T=l 
phases of np-pairing are exclusive at Uu 7^ 0. We believe that this exclusiveness is an artifact 
related to our schematic interaction and/or the assumed self-consistent symmetries. 
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Figure captions 

Fig. 1 Average pairing gaps for a self-conjugate nucleus, N=Z, as a function of x T=0 = 
Gn^/G^p 1 - Left (right) panel shows the results of calculations without (with) particle 
number projection. 

Fig. 2 Calculated additional binding energy, E(x T=0 ) — E(x T=0 = 1), arising from the 
presence of T=0 pairing for a number of Cr-isotopes in the vicinity of the T^=0 line. Different 
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curves denote results for: x T=0 =l.l (•), 1.2 (A), 1.3 (solid triangles) and 1.4 (*). The solid 
line without symbols denotes the Wigner energy term due to [[12|] and the dotted line (*) 
marks the result of the ETFSI-model ||13|| . 

Fig. 3 The calculated dynamical moments of inertia for 46 Cr and 48 Cr. Open circles 
correspond to the case of pure T=l, T^=±l pairing and the sharp peak in is due to 

the breaking of the fr/2 pairs. The curve marked by solid dots indicate calculations with 
undercritical T=0 np-pairing strength at hou—0. The sudden rise of the moment of inertia 
corresponds to the critical frequency, where the T=0 pairing correlations switch on. Note 
the entirely different behaviour of the moment of inertia for the two cases. The dotted line 
(*) denotes the calculations without pairing. 
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